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Abstract — In a multi-agent system, unconditional (multiple) 
consensus is the property of reaching to (multiple) consensus 
irrespective of the instant and values at which states are initial- 
ized. For linear algorithms, occurrence of unconditional (multi- 
ple) consensus turns out to be equivalent to (class-)ergodicity of 
the transition chain (An). For a wide class of chains, chains with 
so-called balanced asymmetry property, necessary and sufficient 
conditions for ergodicity and class-ergodicity are derived. The 
results are employed to analyze the limiting behavior of agents' 
states in the JLM model, the Krause model, and the Cucker- 
Smale model. In particular, unconditional single or multiple 
consensus occurs in all three models. Moreover, a necessary 
and sufficient condition for unconditional consensus in the JLM 
model and a sufficient condition for consensus in the Cucker- 
Smale model are obtained. 

Index Terms — Unconditional single or multiple consensus, 
ergodicity or class-ergodicity, JLM model, Cucker-Smale model, 
multi-agent systems. 

I. Introduction 

Linear consensus algorithms in multi-agent systems have 
gained an increasing attention in various fields of research. 
A general linear consensus protocol is described by the 
following equation: 



X{n+l) = AnX{n), n>0, 



(1) 



where X{n) stands for the vector of states and (A„), the 
transition chain, is a sequence of row stochastic matrices, 
called transition matrices. A matrix is called row stochastic 
if each of its rows sums to 1. The transition matrices may 
be known a priory (exogenous models) or may depend on 
agents' states (endogenous models). 

A simple exogenous model, when the transition chain 
consists of identical row stochastic matrices, was first in- 
troduced by DeGroot in [l]. Later, the same consensus 
problem, but with time dependent transition matrices, was 
considered in 121. The authors related the consensus problem 
to ergodicity of chain (A„) and studied backward products 
of row stochastic matrices to analyze ergodicity of the chain. 
The results of |2| were improved in ||3], ||5], as weaker 
sufficient conditions for consensus were derived. Sufficient 
conditions obtained in |i3l, H, IS) can be expressed briefly in 



two important conditions; non vanishing interaction rates and 
repeated connectivity of communication graph. Alternatively, 
Viscek et al. in |6| modeled a system of agents moving on 
a plane with the same speed but different directions. In the 
Viscek model, the parameter to be updated is the heading of 
agents. Later, by linearizing the Viscek model, Jadbabaie et 
al. |7| obtained sufficient conditions for consensus to occur. 
In the JLM model, each agent updates its state by taking 
the arithmetic mean of the current states of its neighbors 
and itself. It is worth mentioning that the JLM model is 
an exogenous model, since for each agent i, the set of i's 
neighbors at any time is pre-defined. Many authors have 
considered the JLM model and used various techniques to 
weaken the sufficient conditions for consensus as much as 
possible |8|, |9|, [TOl, Ol, [T?]. Together with the notion of 
consensus, existence of individual limits for agents' states, 
i.e., multiple consensus, has been studied in some articles 
such as {131 . lfT4ll . ifTSl ). Moreover, necessary conditions for 
consensus were investigated by Touri and Nedic 1 16J . LITJ . 
ifTSll . In the recent work |fT9l . the authors obtained necessary 
and sufficient conditions for consensus and multiple consen- 
sus for a broad class of consensus algorithms. 

On the other hand, the Umiting behavior of agents in en- 
dogenous model has been widely studied in the literature. The 
Krause model [20] is a well-known example of endogenous 
averaging algorithms. As in the JLM model, each agents 
modifies its states to the arithmetic mean of its neighbors' 
states and its own state. However, unlike the JLM model, 
neighbors of an agent is not pre-defined. More precisely, 
two agents are called neighbors at time n > if their states 
differ by a number less than or equal to some pre-specified 
constant R. In ll20l . the limiting behavior of agents' states 
was studied and it was proved that agents eventually separate 
to disjoint clusters with inter-cluster distances greater than R. 
The Krause model was further studied in lISTI . ll22l . and the 
results were improved. 

Another interesting endogenous linear model was intro- 
duced in ||23]| . known as the Cucker-Smale model or the C-S 
model, in which birds are considered as agents seeking to 
fly with the same velocity in the space. In the C-S model, 



velocities are updated according to an averaging algorithm, 
where birds' interaction rates are functions of distance be- 
tween birds. In [231, sufficient condition for converging to 
the same velocity were derived. 

Our aim is to investigate the limiting behavior of agents in 
various linear averaging models. This paper is organized as 
follows. Section II explicitly describes important notions and 
terminology required for the analysis of limiting behavior. 
The main theorems about consensus and multiple consensus 
are presented in Section III. By considering well-known 
models in Section IV, applications of the main results are 
illustrated. Concluding remarks in Section V close the paper. 

A. Notations 

Through this article, 

• S is the set of agents and s — \S\ is the number of 
agents. 

> n stands for discrete time index. 

. X{n) = [Xi{n) ••• Xs{n)]'^, n > 0, is the state 
vector. 

• Unless otherwise specified, a chain will designate a 
sequence of row stochastic matrices. 

• An, n > is the matrix of interaction rates aij{n), 
l<i,j <s. 

II. Notions and Terminology 

A. Ergodicity and Class-Ergodicity 

Definition 1: Consider a multi-agent system with updates 
equation described in Eq. (HJ. Unconditional consensus in 
system Eq. ([T]) is defined as occurrence of consensus, no 
matter at what instant or values states are initialized. 
We now recall the definition of ergodicity from [16 1. Let 
{An) be a chain of row stochastic matrices. For fc > ^ > 0, 
let A{n, k) = A„_i A„_2 ■■■Ak- Now, 

Definition 2: A chain (A„) of row stochastic matrices is 
called ergodic if for each fc > 0, lim„__j.oo A{n, k) exists and 
is equal to a matrix with identical rows. 
One can show that occurrence of unconditional consensus 
in system Eq. ^ is equivalent to ergodicity of chain (An). 
This describes the relation of unconditional consensus and 
ergodicity. In the following, we define multiple consensus, 
which is another important notion in multi-agent systems. 

Definition 3: Consider the system with dynamics de- 
scribed by Eq. ([T]i. By unconditional multiple consensus in 
system Eq. ([U, we mean that for every i, I < i < s, 
liuin^oo Xi{n) exists, no matter at what instant or values 
states are initialized. 

To formulate multiple consensus as a property of chains of 
row stochastic matrices, we shall introduce class-ergodicity: 
Definition 4: A chain (A„) of row stochastic matrices is 
said to be class-ergodic if for every fc > 0, lim„_j.oo A{n, fc) 
exists and can be relabeled as a block diagonal matrix with 
constant rows. 

Clearly, Class-ergodic of chain {An) implies unconditional 
multiple consensus in system Eq. ([Til. The converse is true 
also Iil9i . Therefore, occurrence of unconditional multiple 



consensus in a system with dynamics described by Eq. ([TJ 
is equivalent to class-ergodicity of chain {An). 

We shall investigate ergodicity and class-ergodicity of 
chains. We now provide essential notions that are employed 
to obtain our main theorems. 

B. li-approximation U7\l 

Definition 5: Chain {An) is said to be -approximation 
of chain (i?„) if 

oc 

^||A„-i3„|| <oo, (2) 

n=0 

with respect to the max norm which is equal to the maximum 
of the absolute values of the matrix entries. 

It is not difficult to show that Zi-approximation is an equiv- 
alence relation in the set of chains of stochastic matrices. 
Importance of the -approximation notion comes from the 
following lemma that is a result of Lemma 1, stated and 
proved in |17|. 

Lemma 1: Let {An) be /i-approximation of chain 
Then, {An) is class-ergodic if and only if (i?„) is. 

C. Infinite Flow and Absolute Infinite Flow 

According to fT6l, chain {An) has the infinite flow prop- 
erty if for every T C S — {l,...,s}, excluding the empty 
set and S the following holds: 

oo 

KjH +aj»(")) = oo- (3) 

n=0 ieT.j^T 

The absolute infinite flow property of a chain is defined as 
follows [18|. 

Definition 6: A chain {An) of row stochastic matrices has 
the absolute infinite flow property if 

where the T(0), . . . is an arbitrary sequence of subsets 
of S with the same cardinality, and Ti denotes the comple- 
ment of Ti in S. 

Clearly, the infinite flow property is implied by the absolute 
infinite flow property. In ||T6l , it has been shown that the 
infinite flow property is necessary for a chain to be ergodic 
[,16J. Following [16j . the authors later showed that the 
absolute infinite flow property is also necessary for ergodicity 
In addition, they proved necessity and sufficiency of 
the absolute infinite property in case of chains of doubly 
stochastic matrices. 

D. Balanced Asymmetry 

We now define another property of chain of row stochastic 
matrices, called balanced asymmetry, which plays a key role 
in our analysis. 

Definition 7: Let {An) be a chain of row stochastic ma- 
trices. Chain {An) is said to be balanced asymmetric if there 
exists a finite M > 1 such that for any two non empty subsets 



Si and S2 of S = {1, . . . , s} with the same cardinality, we 
have 

Remark 1: For those chains that are Zi-approximation of 
balanced asymmetric chains, the infinite flow property is 
equivalent to: 

00 

^^^a,j(n)=oo (6) 

for any subset T of S* as in Eq. (|3]l. Similarly, the absolute 
infinite flow property is equivalent to: 

00 

^ ^ ^ ay(n) = 00 (7) 

n=OigT(n+l)iGT(n) 

for any sequence T{n) of subsets of T as in Eq. This 
can be easily seen by combining Eqs. (|4]i and (|5]l. 

E. Self-Confidence, Type-Symmetry, and 
subsymmetry 

Definition 8: A chain (A„) of stochastic matrices is said 
to have the self-confidence property, or self-confident, if there 
exists ^ > such that for every i, an > S. 
It is worth noting that if a chain is self-confident, then the 
absolute infinite flow property is equivalent to the infinite 
flow property. We also note that for self-confident chains, 
the balanced asymmetry property becomes equivalent to a 
simpler property that is defined in the following. 

Definition 9: ifTSl A chain of row stochastic matrices is 
type-symmetric if there exists a finite M > such that for 
every rc5 = {l,...,s}, excluding the empty set and S, 

In the following, according to [14], we define a special 
case of type-symmetry. 

Definition 10: A chain (A„) of row stochastic matrices is 
called sub-symmetric if there exists a finite AI > such that 
for every two agents i and j, we have 

atjin) < Maj,{n), Vn > 0. (9) 

F. Strong Interaction Digraph of A Chain 

The strong interaction digraph of a chain is of importance 
in class-ergodicity analysis. To evaluate the total influence of 
an agent on another agent, we define a function of ordered 
pairs of agents, as follows. 

Definition 11: Let [An) be the chain representing interac- 
tion rates between s agents. We define function int : SxS ^ 
R^" by 

00 

mt(«, j) = ^ ay (n). (10) 

n=0 

An ordered pair (i,j) of agents is said to be highly interactive, 
if 

00 

mt(i,j) = ^ aij(n) = 00. (11) 

Tl = 



Note that from being highly interactive, one cannot 

conclude that {j, i) is highly interactive. 

For a chain (A„), strong interaction digraph Ga is formed 
as follows. Considering (A„) as interactions between s 
agents, Ga contains s nodes, each representing an agent. 
An edge is drawn from node i to node j if and only if the 
ordered pair is highly interactive. 

One defines a relation R between the s nodes of Ga 
(agents) as follows. For every two nodes i and j that are not 
necessarily distinct, iRj if and only if there is a directed 
path from i to j. 

The following lemma, as stated in lil9J . shows how strong 
interaction digraph is related to balanced asymmetry. 

Lemma 2: Assume that {An) is a balanced asymmetric 
chain, and Ga is its strong interaction digraph. Then, R is 
an equivalence relation 

Note 1: The equivalence relation R partitions agents into 
equivalence classes, herein called islands. 

III. Main Theorems 

In this section, we state our main results in two parts: 
ergodicity results and class-ergodicity results. 

A. Ergodicity Results 

The following theorem is established in ||T9l . 

Theorem 1: If chain (A„) is -approximation of a bal- 
anced asymmetric chain, then (A„) is ergodic if and only if 
it has the absolute infinite flow property. 

B. Class-Ergodicity Results 

From [19 1, we have the following theorem on class- 
ergodicity of chains. 

Theorem 2: Let chain (A„) be Zi-approximation of a bal- 
anced asymmetric chain. Then, chain {An) is class-ergodic if 
and only if the absolute infinite flow property holds over each 
island of the strong interaction digraph induced by {An). 

In the following, by employing Theorem |2] we prove a 
sufficient condition for a chain to be class-ergodic. 

Theorem 3: If chain {An) is /i-approximation of a self- 
confident and type-symmetric chain, it is also class-ergodic. 

Proof: From Lemma [T] to prove class-ergodicity of 
{An), we can assume that {An) is self-confident and type- 
symmetric. These two properties of {An) imply that (A„) 
is balanced asymmetric. Therefore, according to Theorem 
it suffices to show that the absolute infinite flow property 
holds over each island of the strong interaction digraph 
Ga- Let / be an arbitrary island and T(0),T(1), ... be an 
arbitrary sequence of subsets of / with the same cardinality. 
Calling the equivalence of Eqs. (|7]i and (|4]l for chains that 
are /i-approximation of balanced asymmetric chains, in the 
following, we show that Eq. ^ holds by considering the 
following two cases: 

Case I. The sequence T(0), T(l), . . . becomes constant 
after a finite time, i.e., there exist T c I and > such 



that T{i) — T for every n > iV. In this case. 



^ ^ ^ a.,H> 5]^5]a,,(n) (12) 



n=0 ieT(n+l) jeT{n) 



n=N ieTjeT 



Clearly, there exist two agents p E T and q <E T such that 
pair (p, q) is highly interactive. Otherwise, island / would 
not be strongly connected. Noting that the RHS of ( fT^ is 
not less than 



^ apq{n) 



(13) 



n=N 



which diverges, since (p, q) is highly interactive, the result if 
proved. 

Case 11. The sequence T(0),T(1), . . . does not converge, 
i.e., there exists a time subsequence no,ni,... such that 
T(nfc) 7^ T{nk + 1) for every fc = 0, 1, . . .. Clearly, 

Sn=0 Sief (n+1) X]jeT(n) '^ij (") 

— Sfc=0 SieT(nfc + l) X]jeT(nfc) 0-ij{nk) 

Since T{n}^) ^ r(nfc + 1) and the two subsets are of the 
same cardinality, there exists an agent that belongs to both 
T{nk + 1) and Tijik). Hence, due to self-confidence of chain 
{An), we have 



(14) 



iGT(nfc + l)ieT(nfc) 



(15) 



Therefore, the RHS of Eq. (fT4l i is not less than X^fc^o ^ 
diverges. This proves the result again in this case. ■ 

IV. Relationship to known models 

In this section, we apply our main theorems to chains 
corresponding to different types of models and consensus 
algorithms found in the literature in order to analyze when 
their transition chains become ergodic or class-ergodic. 

A. JLM Model 

The parameter considered in Q is the heading of each 
agent. If we write 9i (n) as the heading of agent i at moment 
n, the model describing evolution of headings is 



^in + l) 



1 



l + d,{n) 



mn) 



jeDi(n) 



(16) 



where Di{n) and di{n) denote respectively the set and the 
number of neighbors of the agent i at time n. It is also 
assumed that for each n > 0, if j G Di{n), then i e Dj{n) 
too (undirected graph). 

In Q, the authors proved that a sufficient condition for 
consensus to occur is existence of an infinite sequence of 
contiguous, nonempty, bounded, time-intervals [n^, rii+i), 
« = 0, 1, . . ., such that across each such interval, the s agents 
are linked together 

In the following, we wish to apply Theorems [T] and |3] to 
the transition chain of the JLM model. To take advantage of 
Theorem |3] we show that in the JLM model, the transition 
chain is both self-confident and type-symmetric. Note that 



(17) 



This proves the self-confidence of the chain. To prove the 
type-symmetry, it suffices to show that the chain is sub- 
symmetric. If at time n > 0, we have j ^ Di{n), then 
i ^ Dj{n) either. Therefore aij{n) = aji{n) = (con- 
sistently with the subsymmetry requirement). If j e Di{n) 
then i e Dj{n) also. In this case, it is easy to see that aij{n) 
and aji{n) both lie in the interval [1/2, 1/s]. Therefore, the 
subsymmetry condition holds by setting AI = s/2. Thus, 
from Theorem |3] we conclude that the chain is class-ergodic. 
In other words, in the JLM model, unconditional multiple 
consensus occurs without any additional assumption. 

We also note that the chain is balanced asymmetric as 
well, since self-confidence and type-symmetry imply bal- 
anced asymmetry. Thus, from Theorem [T] we obtain that the 
absolute infinite flow property is equivalent to the ergodicity 
of the chain. On the other hand, since the chain is self- 
confident, the absolute infinite flow property is equivalent 
to the infinite flow property. Hence, in the JLM model, 
the infinite flow property is necessary and sufficient for 
ergodicity of the transition chain. 

Graph interpretation of the infinite flow property in the 
JLM model is as follows. Due to the subsymmetry property 
of the JLM model, if a pair is highly interactive, 

then (j, i) is highly interactive also. Therefore, in this case, 
the strong interactions digraph can be considered as an 
undirected graph. The infinite flow property is now equivalent 
to connectivity of the strong interactions graph. 

Another equivalent condition to the infinite flow property, 
that is more similar to the condition derived in [7], is exis- 
tence of an infinite sequence of contiguous and non empty 
time-intervals [n^, i > 0, with the property that across 

each such interval, the s agents are linked together Note that 
the boundedness of the time-intervals is not required unlike 
in Q. More importantly, the condition derived here is not 
only sufficient, but also necessary for ergodicity of the chain 
(unconditional occurrence of consensus in the model). On the 
other hand, unlike |7|, without extra conditions, no statement 
can be made about the speed of convergence to consensus. 

B. Models with Finite Range Interactions 

The Krause model |20| is an example of endogenous 
models with finite range interactions. These models are 
special cases of first order models in which interaction rates 
depend directly on states. In these models, agent i receives 
information from agent j if and only if the distance between 
the two agents is less than some pre-specified level Ri, which 
is in general different for distinct agents. In the following, we 
define the interactions rates between agents. For every agent 
i, we set a decaying function fi : R-^ — i?-" that vanishes 
at Ri, and define 

We consider a particular case of models defined above. 
Assume that agents have the same range of connectivity, i.e., 
Ri = R for every i, and use identical decaying functions, 
i.e., fi = f. In the Krause model, f{x) = 1 for < a; < _R 
and f{x) — elsewhere. 



It can be proved that in this case, the transition chain is 
self-confident with S = 1/s. It can also be shown as follows 
that the transition chain is sub-symmetric. If at time n > 0, 
agents i and j do not communicate, then aij (n) = aji (n) = 
0. If the two agents communicate, then f(\xi{n) ~Xj{n)\) > 
0. Using Eq. ( fTsl l. we have 

^'^jin) ^ J2l=ifi\Xjin) - Xk{n)\) ^^^^ 
^A^) Efc=i fi\X^{n) - Xk{n)\) 

Noting that / is non increasing and f{\Xi{n) — Xi{n)\) = 
fi\Xj{n) - Xj{n)\) = /(O), we conclude that the RHS 
of Eq. ( fT9] l lies in interval [l/s,s]. Hence, subsymmetry 
is established by setting M = s. The chain being both 
self-confident and sub-symmetric, it is also type-symmetric. 
Thus, according to Theorem|3] the chain is class-ergodic, i.e., 
unconditional multiple consensus occurs. 

C. The C-S model 

The C-S (Cucker-Smale) model ["231 is an example of en- 
dogenous consensus models with interaction rates remaining 
strictly positive. We apply our results to a generalized version 
of the C-S model |23 | that describes evolution of positions 
X,;'s and velocities V^'s in a bird flock, in a three dimensional 
Euclidian space: 

X,{n + l) = X,{n) + hV^{n), 

V,{n + 1) =V,{n)+Y.,^J-{\\X,{n)-X,{n)\\) (20) 

{y,{n)-Vdn)), 

where / : R-" — > is a non increasing function. Note 
that in this model, the limiting behavior of velocities is of 
interest. We have 



a,,{n)= S{\\X,{n)^X,{n)\\), \/i ^ j 



(21) 



and 



au{n)^l-yj{\\X,{n)-Xj{n)\\), Vz (22) 



Clearly, the transition chain in this algorithm is symmetric. 
To enforce self-confident, one may require an additional 
assumption, such as, f{y) < 1/s for any y > 0. By this 
assumption, we have aii{n) > 1/s for every i — l,...,s 
and n > 0. Therefore, the chain becomes self-confident. 
The combination of the self-confidence and type-symmetry 
of the chain allows an application of Theorem |3] to yield the 
following result. 

Theorem 4: Consider the algorithm described by Eq. (|20t . 
If f{y) < l/'^ for any y > 0, then the transition chain 
is class-ergodic, and consequently, unconditional multiple 
consensus occurs. 

To state the consensus result for the generalized C-S model, 
we define parameters Mx and My calculated from initial 
positions and velocities: 



Mx 



max{||X,(0) - Xj{0)\\\l < i < j < s} (23) 



max{||F,(0)-F,(0)|||l<z<j<s}. 



(24) 



Theorem 5: For the multi-agent system with dynamics 

described by Eq. (l20l i. assume that f{y) has the following 
property: 

f{y) < 1/s, Vs > (25) 



Assume also that 



(26) 



Then, all agents' velocities converge to a common value. 
Moreover, there exists a non negative number R such that 
for every 1 < i,j < s, 



\X,in) ~ Xj{n)\\ < i?, Vn > 



(27) 



Unlike the models described previously. Theorem |5] is not an 
immediate result of Theorems [T] and |3] However, to prove 
Theorem |5] we employ a technique similar to that used in 
the proof of Theorem [T] in [19|. 

Proof: For every z = 1, . . . , s, let V^^ (n), Vi2 (n), V^a (n) 
be components of Vi{n), i.e., 

V,in) = [Vain) Va{n) V,^{n)]^ (28) 

It is straight forward to verify that V^^'s {r ~ 1, 2, 3) satisfy 
for straightforwardly identifiable coefficients aij, the same 
update equation as V^'s do, i.e.. 



Vr{n+1)- Vr{n) - 2^ a,j{n){Vjr{n) - VM) (29) 
One can rewrite Eq. (|29] l as the following: 

I VJn^ 

(30) 



Vir{n + 1) = [l- '^j^i ^tr{n) 



Let us define Zir{n) : R-" — >-R, l<j<s, r=:l,2,3 from 
Vir{nys, as follows. At every time n > 0, Zir{n) is equal to 
the jth least number among Vir{Ti), . . . , Vsr{n). 

Note that the coefficients in the RHS of Eq. ( |30] | are all 
positive and add up to 1. This means that Virin + 1) is a 
convex combination of the values Vir{n), . . . , Vsr{n). Thus, 
for every r, 1 < r < 3, interval [zir{n), Zsr{n)], which 
is the smallest interval containing all the values Vir{nys, 
shrinks during time. Particularly, this shows that Vir{n), and 
consequently Zir{n), are bounded for every i = 1, . . . , s and 
r = 1,2,3. For consensus to occur, we require that interval 
[zir{n), Zsr{n)] converges to a point for every r, 1 < r < 3. 
In the following, we investigate how these intervals shrink 
with time. Let us define: 

3 



We know that 

X,{n) - X,{n) = X,{Q) ~ X,{Q) 

Thus, 

\\X,{n) - X,{n)\\ < \\XM-X,m 



(31) 



(32) 



(33) 



However, 

\\V^{m) - Vj{m)\\ < J2l=i \V^r{m) - V,rim)\ 

< J2l=ii^sr{ni) - zir{m)) (34) 
= z{m) 

Eqs. (I23]l, (l33]l and (O imply 

n-l 

\\X,{n) - Xj{n)\\ <M, + hJ2 (35) 

rn—O 

From / being non increasing, we obtain 

f{\\X,{n) - X,{n)\\) > f{M, + ^ z{m)) (36) 

Note that the RHS of Eq. ( [36] l above is independent of i, j. 
Defining 



(n) = + ^ z(to), 



5 

Eq. ( |36] l impHes 

atj{n) > figin)), Vn > 



(37) 



(38) 



We know recall Eq. (|30] l. It is straight forward to verify 
that all the coefficients in the RHS of Eq. ( [30l l lie between 
f{g{n)) and 1 — (s — l)f{g{n)). Thus, as the sum of 
coefficients is 1, to find a lower bound for the value of 
Vir{n + 1), we put higher weights on lower valued Vj,/&, 
in particular, 1 — (s — l)/(g(n)) on the least one, which is 
zir{n) and f{g{n)) on the rest of them. Hence, we conclude 



V,rin + 1) > (l-(.s-l)/(.g(n)))zi,(n) 
+ J2j=2fi9in))zjr{n) 

As a result, 

zir{n) > (l- {s - l)f{g{n))jzir{n) 



(39) 



(40) 



Using the opposite process to build an upper bound for the 
value of Vir{n + 1), we obtain: 



Zsr{n) < (^1 - (s - l)f{g{n))jzsrin) 



(41) 



(42) 



Subtracting Eq. (|40] | from Eq. dTTJ implies 

Zsr{n + 1) - Zir{n + 1) < 

(l - sf{g{n))j (^Zsr{n) - zir{n) 

By adding up Eq. ( l42b for r = 1, 2, 3, we obtain 

z(n + l) < (l-.s/(g(n)))z(n) (43) 
or equivalently, 

z{n + 1) - 2:(n) < -sf{g{n))z{n) (44) 



We now note that hz{n) — g{n + 1) — g{n). Since / is non 
increasing and g{n + I) — g{n) ~ hz{n) > 0, we have 



f{g{n))z{n) = f{g{n))Sill±n^ 



(45) 



1 f9(«+l) 

Eqs. dllll and (|45]l imply 

z(n + l) -z(n) < — / f{y)dy (46) 

The above equation holds for every n > 0. If we sub- 
stitute variable n in Eq. ( l46l l with n' and sum it up for 
7i' = 1, . . . , n — 1. We obtain 



z(n) - z(0) < 



9(") 



I{v)dy 



(47) 



9(0) 



Recalling the definition of g{n) we conclude 

z(n) - z(0) < — / f{y)dy (48) 

If consensus does not occur, then Em=o '^(") diverges. Thus 
by taking n to infinity and noting that lim„_^tx) z{n) exists as 
z(n) is non increasing and non negative, if consensus were 
to fail, Eq. ( l48T l would be modified as 

lim z(n) - z(0) < — / f{y)dy (49) 

On the other hand, according to our assumption Eq. (l26i . we 
know that 

^(o) = E'=i(^-(o)-zi.(o)) 

= E'=iniax,,,{(F,,(0)-y,.(0)} 
< Er=imax,j{(||t/,(0) - y,(0)||} = 3M„ 

Eqs. (|49] l and (fSOl l together imply that 

lim z(n) < (51) 

which is a contradiction, since z(n) is non negative. Hence, 
consensus must occur and moreover, EJ^i z{n) converges. 
Recalling Eq. ( [35l ) we conclude that ||Xi(n) — Xj{n)\\ is 
bounded for every i^j, i.e., there is i? > such that 

\\X,{n) - Xj{n)\\ < R, Vi, j, Vn > (52) 



(50) 



Applying Theorem |5] to the C-S model with 



(53) 



we obtain the following. 

Corollary 1: Let the dynamics of a multi agent system be 
described by Eq. (|20] | with / defined in Eq. (l53b . Assume that 
K/a"^^ < 1/s. Then, under any of the following conditions, 
agents velocities converge to a common value: 

1) /?< 1/2 

2) ^ > 1/2 and 

sK 



M„ < 



3ft(2/3-l)(M, + a)2/3-i 



(54) 



V. Concluding Remarks 

In this paper, we analyzed the Umiting behavior of the JLM 
model, the Cucker-Smale model, and the Krause model, using 
our general theorems about consensus and multiple consen- 
sus. In Theorem [T] we presented a necessary and sufficient 
condition for chains that are ^i-approximation of balanced 
asymmetric chains to be ergodic. Noting that the transition 
chain in the JLM model is balanced asymmetric, we found a 
necessary and sufficient condition for unconditional consen- 
sus to occur in the JLM model. Theorem |2] provides a nec- 
essary and sufficient condition for class-ergodicity of those 
chains that are ^i-approximation of balanced asymmetric 
chains. Theorem [3] is a non trivial especial case of Theorem 
121 that contains a sufficient condition for class-ergodicity. 
Theorems |2] and |3] led us to prove unconditional multiple 
consensus in various known models: the JLM model, finite 
range interactions models, the Krause model for instance, 
and the generalized version of the C-S model. Furthermore, 
we obtained a sufficient condition for consensus in the C-S 
model by employing elementary methods. 

In future work, we wish to investigate the rate of con- 
vergence when consensus or multiple consensus occurs in 
a system. In particular, we are interested in necessary and 
sufficient conditions for exponentially fast convergence. We 
also wish to extend our results to random networks as well 
as systems with large population of agents. 
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